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ABSTRACT 

The physical origin of radio emission in Radio Quiet Active Galactic Nuclei (RQ AGN) 
remains unclear, whether it is a downscaled version of the relativistic jets typical 
of Radio Loud (RL) AGN, or whether it originates from the accretion disk. The 
correlation between 5 GHz and X-ray luminosities of RQ AGN, which follows = 
lQ~^Lx observed also in stellar coronae, suggests an association of both X-ray and 
radio sources with the accretion disk corona. Observing RQ AGN at higher (mm-wave) 
frequencies, where synchrotron self absorption is diminished, and smaller regions can 
be probed, is key to exploring this association. Eight RQ AGN, selected based on their 
high X-ray brightness and variability, were observed at 95 GHz with the GARMA 
and ATCA telescopes. All targets were detected at the 1 — 10 mJy level. Emission 
excess at 95 GHz of up to x7 is found with respect to archival low-frequency steep 
spectra, suggesting a compact, optically-thick core superimposed on the more extended 
structures that dominate at low frequencies. Though unresolved, the 95 GHz fluxes 
imply optically thick source sizes of 10“"* — 10“^ pc, or ^ 10 — 1000 gravitational 
radii. The present sources lie tightly along an La (95 GHz) = 10“"^Lx (2-10 keV) 
correlation, analogous to that of stellar coronae and RQ AGN at 5 GHz, while RL 
AGN are shown to have higher Ln/Lx ratios. The present observations argue that 
simultaneous mm-wave and X-ray monitoring of RQ AGN features a promising method 
for understanding accretion disk coronal emission. 
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1 INTRODUCTION 


Radio loud (RL) active galactic nuclei (AGN) are 
known for their relativistic jets that often extend be¬ 
yond the nucleus. Radio e mission from Rad io Quiet 

(RQ) AGN (as defined by iKellermann et all llQsJ l is 
lower by several orders of magnitude than that of 
RL AGN. RQ AGN show, on average, smaller radio 
structures some of which are often unresolved down 
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scaled version of the RL, collimated jets llBarvainis et al.l 
1 19961 : iGallimore et al.ir2006l l , or coronal emi ssion from mag¬ 
netic activity above the accretion disk (IField Roeersl 
I1993II , remains to be resolved. In some sources, a com¬ 
bination of several spectral components may be present 
aroun d the nucleus dBarvainis et al.l 1 19961 : IGallimore et al.l 
I2OO4II . For the sake of clarity of discussion, in this paper we 
distinguish a jet that is both highly relativistic and well- 
collimated, from coronal emission that arises from hot gas, 
and perhaps an outflow that is neither relativistic nor well- 
collimated. 


A connection between the radio and X-ray emission 
in RQ AGN is suggested by the correlation of the ra¬ 
dio luminosity La at 5 GHz and the X-ray luminosity 


I2OO6 ; Panessa et al.ll2007ll. Laoi_&BehaiJ ( 2008f ) used the 
PG quasar sample Jschmidt &: Green 1983ll to demonstrate 
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that not only are L[i{= uLn at 5 GHz) and Lx (0.2 - 20 keV) 
correlated over a large range of AGN luminosity, but that the 
correlation follows the well establ ished relation for cor onally 
active cool stars Lr/Lx ~ 10“® llGiidel fc Benj[l993l ). This 
suggests that radio emission from RQ AGN is due to mag¬ 
netic coronal activity, akin to that of stellar coronae. Since 
the LrILx ~ 10“® relation is accepted as a manifestation 
of coronal heating by energetic electrons following magnetic 
reconnection in cool-stars, which subsequent ly gives rise to 
X-ray emission, the correlation presented in iLaor fc Behatl 
ll2008ll over 20 orders of magnitude in luminosity raises the 
possibility that radio emission in RQ AGN may also be re¬ 
lated to coronal, magnetic activity. 

There are, however, several physical differences between 
X-ray and radio emission of stellar coronae and AGN. X- 
ray spectra of stellar coronae are thermal with T ~ 10® K, 
while the non-thermal (Comptonized) spectra of AGN imply 
T ~ 10® K. Moreover, while X-rays from AGN vary dramat¬ 
ically over short time scales, the variability at 5 GHz and 
8.5 GHz is much slower and smaller in amplitude. If radio 
emission from RQ AGN is due to self absorbed synchrotron, 
absorption should decrease with frequency. Consequently, 
observations at higher frequencies should probe smaller re¬ 
gions of the source, and variability time scales of RQ AGN 
should thus also dec rease with fr e quenc y. Very recently, a 
theoretical work by llnoue fc Doil (l2014l . arXiv 1411.2334) 
further highlights the potential of mm observations for de¬ 
tecting coronal magnetic activity on AGN accretion disks. 

As opposed to 5 and 8.5 GHz, higher frequenc y observa¬ 
tions of RQ AGN are scarce. iBarvainis et al.l(ll996ll observed 
RQ AGN up to 15 GHz. About half of the RQ AGN in that 
sample appeared to have flat or inverted spectral compo¬ 
nents, and a hint of variability, strongly suggesting o paque 
synchrotron emission from the nucleus. IPark et al.l (l2013l l 
reach similar conclusions based on observations at 22 and 
45 GHz. Above 300 GHz (sub-mm) there is a steep rise of 
the spect rum due to the low frequency tail of therm al dust 
emission (IBarvainis et al. I I 1992 I : ikughes et al.l Il993l l , which 
leaves a largely unexplored band between 50 - 300 GHz, 
crudely referred to here a s the mm waveband. N otable ex¬ 
ceptions are the works of iDoi et ahl ll2005l . l201l[l , who ob¬ 
served low luminosity AGN and early type galaxies. That 
sample includes both RQ and RL sources, some of which fea¬ 
ture flat or even inverted spectra at high frequencies, indicat ¬ 
ing genuine optically thi ck cor e emission. lOw en et all lll98lf) , 
lOwen fc Puschelj lll982ll , and I Sadler et al.1 ( 2008l ~l studied X- 
ray selected quasars, but mostly RL AGN brighter than 1 Jy 
at 90 and 95 GHz, much brighter and more luminous than 
any of the present sources. 

Little data are available on radio monitoring of RQ 
AGN on short time scales. Even less are availab le for simul¬ 
taneou s radio and X-ray observations. Recently. I Jones et al.l 
(120111) published their campaign for NGC4051. Despite the 
systematic monitoring, the conclusion of that study was 
ambiguous due to the limited radio variability at 8.4 and 
4.8 GH z. It is also worth noting a similar work bv iBell et al.l 
ll201ll) . but for the RL NGC 7213. 

In this work, we focus on eight RQ AGN that are well 
known for their bright X-ray emission and variability. We ob¬ 
serve them at 95 GHz in order to obtain their high-frequency 
radio properties. One aspect we wish to test is whether they 
are suitable candidates for high-frequency radio (and X-ray) 


monitoring. Ultimately, we believe the simultaneous radio 
and X-ray monitoring of these targets will enable the ex¬ 
ploration of the physical origin of the radio emission in RQ 
AGN. If the aforementioned coronal conjecture is valid, at 
high radio frequencies, variability is expected to approach 
and be closely related to that of the X-rays, as is often the 
case in stellar coronal flares, and manifested, e.g., by the 
Neupert effect llNeupertlll968h . The remainder of the paper 
is organized as follows: In Sec. [5] we describe the observa¬ 
tions of the eight objects. In Sec. |3] we present the results, 
including source flux densities and inferred broad-band ra¬ 
dio spectra. The implications of these spectra are discussed 
in Sec. |4] with a comparison with RL AGN. The conclusions 
and outlook are given in Sec. (5] 


2 SAMPLE AND OBSERVATIONS 

The current sample was chosen based on the high X-ray 
brightness and documented variability of its sources. This 
criterion selects high accretion AGN among the RQ ones. 
Although the flux and morphology of most AGN at 95 GHz 
was essentially unknown, we did make sure to choose radio 
sources that are compact and sufficiently bright (> 1 mjy) 
at lower frequencies. Six of the targets, namely NGC 3783, 
NGG5548, NGG7469, NGC 3227, Ark 564, and Mrk766, 
are highly variabl e in the X-rays accord i ng to the exten¬ 
sive monitoring of iMarkowitz fc EdelsonI (l2(304h . MR 2251- 
178 is a bright X-ray quasar that has been observed of¬ 
ten in the X-rays. The X-ray light curve of MR 2251-178 
closely correl ates with the optical light curves over a period 
of 2.5 years (I Arevalo et al.1 [20081) . MR 2251-178 as well as 
NGC 3783 show e vidence for variable emission in the near-IR 
llLira et al.ll201ll ). Finally, NGC 5506 is the brightest radio 
source in the sample. In the X-rays, it re veals X-ray variabil¬ 
ity on time scales from days to years llUttlev fc McHardvl 
I2OO5I: ISoldi et al.ll201lD . 

We observed continuum emission in the 3 mm (95 GHz) 
window from all eight targets. Three targets were observed 
with ATGA (the Australia Telescope Compact Array), and 
hve targets were observed with CARMA (the Combined Ar¬ 
ray for Research in Millimeter-wave Astronomy). The obser¬ 
vation log from both telescopes is given in Table [T] 

ATCA is an array of six 22-m antennas at the Paul Wild 
Observatory (Australia). The ATCA data were taken in 
September 30, 2013 and in March 29, 2014. We used ATCA 
in an extended array configurations (H214D and H168D) to 
provide a good uv coverage with long baselines (up to 4500 
m) and reach an angular resolution of 2"at 95 GHz. 

CARMA is a 15 element interferometer consisting of 
nine 6.1 meter antennas and six 10.4 meter antennas, lo¬ 
cated in California (USA). The CARMA observations were 
performed in C- or D-array configuration, reaching an an¬ 
gular resolution of 2'.'2 and 5”5, resp ectively. _ 

The MIRIAD software package (ISault et al.lll99^ was 
used for reduction of the visibility data, including flagging 
data affected by shadowed antennas, poor weather or an¬ 
tenna malfunctions. Absolute flux, phase, and pass band 
calibrators are listed in Table [T] After obtaining the phase 
and amplitude solution, we apply them to the target source 
using standard procedures. Eventually, we invert and clean 
with natural weight the visibility data to obtain the map. 
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Table 1. Observation log 


Object 

Telescope 

Array 

Resolution 

(") 

Date 

Phase 

calibrator 

Pass band 

calibrator 

Flux 

calibrator 

MR 2251-178 

ATCA 

H214D 

2.0 

30-09-2013 

2243-123 

3C446 

Uranus 

NGC 3783 

ATCA 

H168D 

2.0 

29-03-2014 

1144-379 

3C279 

Mars 

NGC 5506 

ATCA 

H168D 

2.0 

29-03-2014 

1406-076 

3C279 

Mars 

NGC 7469 

CARMA 

C 

2.2 

11-11-2013 

3C454.3 

1927-1-739 

mwc349 

ARK 564 

CARMA 

C 

2.2 

16-11-2013 

2236-1-284 

1927-1-739 

mwc349 

NGC 3227 

CARMA 

D 

5.5 

03-01-2014 

0956-1-252 

3C273 

Mars 

MRK 766 

CARMA 

D 

5.5 

19-01-2014 

1159-1-292 

3C273 

Mars 

NGC 5548 

CARMA 

D 

5.5 

17-02-2014 

1310-1-323 

1635-1-381 

Mars 


The MIRIAD imfit task is used to measure the source flux 
density as well as its uncertainty using a point-source fit. 


3 RESULTS 

3.1 95 GHz Flux Density 

We obtained radio maps for all of the objects observed with 
the CARMA and ATCA telescopes at 95 GHz. The maps are 
presented in Fig. [T] The measured flux densities range be¬ 
tween ~1 to 10 mjy. The sources are essentially unresolved, 
or just resolved at the observed scales (Table [T]). This in¬ 
dicates that the morphology of the radio source, whether a 
narrow relativistic jet, or a broader and slower coronal out¬ 
flow, can not be assessed from imaging. Consequently, the 
essential physical information of the source will currently 
need to come from the broad-band radio spectrum, and later 
from variability. 

The presently measured flux densities at 95 GHz are 
listed in the first column of Table [2] along with their uncer¬ 
tainties, which are generally < 20 %, and much better for 
the brightest sources. The known distances and consequent 
luminosities are also listed (with the same relative uncer¬ 
tainties). These fluxes are likely dominated by core emis¬ 
sion. However, the possible contamination from unresolved 
radio jets cannot be ruled out. The physical size in pc of an 
opaque (self-absorbed) synchrotron source can be estimated 
from the measured flux (e.g., Laor & Behar 2008) to be: 

^pc = (1) 

where izcHz is the observed frequency in GHz, L 30 = 
ergs“^ Hz~^ and B is the magnetic held in Gauss. 
Hence, the size of this optically thick radio-sphere increases 
with radio luminosity, and decreases with frequency. The 
term radio-sphere is used below when referring to this es¬ 
timate of the radio-core size based on its luminosity. The 
measured 95 GHz fluxes of the eight sources, therefore, can 
provide an idea of the sizes of the 95 GHz cores, assuming a 
value B for the magnetic field, although the dependence on 
B^^^ is weak. The resulting sizes of the eight present sources 
are listed in Table[21 assuming B = 1 G, which is of the order 
expected from equipartition (Beq/ 87 r = Lboi/dyrRpcC, where 
I/boi is the bolometric luminosity). A n expression for R^c a t 
equipartition is given in Eq. (22) in iLaor fc Behail (l2008li . 
Ignoring the unknown B, the (fractional) uncertainty of f?pc 
can be assumed to be half that of the flux density (second 
column in Table H, as Rpc oc , namely < 10%. Evi¬ 


dently, the obtained sizes Rpc of the 95 GHz radio-spheres 
are of the order of 10“^^ to 10~® pc. The light crossing time 
Rpc/c can provide an idea of the minimally expected vari¬ 
ability time scales at 95 GHz, namely 0.1 — 1 day. 

The 2-10 keV X-ray luminosities Lx (accurate to a 
few percent) are also listed for comparison in Tabled The 
95 GHz to X-ray luminosity ratios of the present sample are 
rather tightly concentrated around Lh/Lx = 10 ““*, with 
loglLii/Lx) = —4.1 ± 0.25 (Table [21. The chance proba¬ 
bility for such a tight relation is only < 0.00075. How this 
compares with other samples is discussed in Sec. WT\ 

The black hole mass Mbh, Rpc in terms of gravita¬ 
tional radii Vg = GMb,yi/(? , and Lboi/LBdd are listed in 
Table O Evidently, the radio-sphere sizes correspond to 
« 10 — 1000 Tg. The bolometric to Eddington luminosity 
ratio Lboi/LEdd of the present sample also spans more than 
two orders of magnitude from 0.01 - 3.0. These ranges are 
much larger than the uncertainties in MBH(LBdd) and Lboi, 
which ar e dominated by syst ematics, and can reach a fac¬ 
tor of 3 ([Peterson et al.lr2004l l. In the present small sample, 
we confirm that Lboi scales with LgsGHz, whi ch is consis¬ 
tent with low-frequency results (Lboi oc in 


2007) and with radio galaxies (Lboi oc L^, in 


White et al 


Sikora et al 


201311 . 


The factor 10 range in RpRoc Table[2I, and factor 
of 100 range in rg{oi Mbh, Table 0 results in an inverse 
correlation of RpcRg with Mbh. The four low-mass black 
holes with Mbh < lO’^ all have RpcRg ^ 330, and the four 
massive ones with Mbh ^ 3 x 10^ all have Rpc/vg ^ 52. 
Interestingly, RpRrg appears to increase with Lboi/LBdd in 
the sense that the four low Lboi/LBdd(< 0 . 1 ) sources have 
small (weak) radio-spheres of a few lOr^, while the four 
high Lboi/LEdd(> 0.4) sources have large (luminous) radio¬ 
spheres of 300 — 1000 Tg, implying that higher accretion rates 
(Lboi/T/Edd) are associated with mm-wave activity on larger 
relative scales. 


3.2 95 GHz Excess 

We used the literature to obtain the radio flux density of 
the present sources at frequencies below 95 GHz. We found 
measurements from 1.4 GHz up to the highest radio fre¬ 
quency available of 22 GHz for NGC5506 and NGC7469. 
When the core is resolved in the archival data, we only con¬ 
sider that component to probe the emission closest around 
the nucleus. In order to obtain the highest spatial resolu¬ 
tion, we first consider radio detections from the VLA in 
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MR2251-178 NGC3783 NGC5506 



Right ascension (J2000) Right ascension (J2000) Right ascension (J2000) 


NGC7469 ARK564 NGC3227 



23 0316.0 15.8 15.6 15.4 15.2 15.0 22 42 39.8 39.6 39.4 39.2 39.0 38.8 10 2331.5 31.0 30.5 30.0 

Right ascension (J2000) Right ascension (J2000) Right ascension (J2000) 


MRK766 NGC5548 



1218 27.2 27.0 26.8 26.6 26.4 26.2 26.0 25.8 1418 00.4 00.0 59.6 59.2 58.8 

Right ascension (J2000) Right ascension (J2000) 


Figure 1. 95 GHz contour maps from the present observations. Beam size is indicated for each observation. Contour levels of the maps 
are: for MR 2251-178 (-0.4,0.4,0.8,1.2) mjy/beam; for NGC 3783 (-0.8,1.2,1.9,2.3,2.9) mjy/beam; for NGC 5506 (-1.3,2.5,3.8,6.3,7.5) 
mjy/beam; for NGC7469 (-0.5,1,2,3,4) mjy/beam; for ARK 564 (-0.38,0.54,0.9,1.1) mjy/beam; for NGC 3227 (-0.4,1,1.7,3,3.8) 
mjy/beam; for MRK 766 (-0.4,0.6,1.0,1.2,1.8) mjy/beam; for NGC 5548 (-0.6,1,1.4,1.6) mjy/beam. 


Table 2. Luminosity and radio-sphere size properties 


Object 

GHz 

mjy 

Ol 

Mpc 

Lgs GHz 

lO^'^ergs^^Hz”^ 

Rpc 

10-3 pc 

log Lr,(95 GHz) 

erg 

log Lx 

erg 

log{Lii/Lx) 

MR2251-178 

1.6±0.3 

289.1 

0.16±0.03 

0.70 

40.18 

44.40 

-4.22 

NGC 3783 

3.1±0.7 

42.2 

0.0066±0.0015 

0.15 

38.78 

43.15 

-4.35 

NGC 5506 

lOitl.O 

26.8 

0.009±0.001 

0.17 

38.91 

42.74 

-3.83 

NGC 7469 

4.95±0.16 

71.2 

0.030±0.001 

0.32 

39.46 

43.18 

-3.72 

ARK 564 

1.14±0.19 

108.4 

0.016±0.003 

0.23 

39.18 

43.57 

-4.39 

NGC 3227 

4.1±0.24 

16.7 

0.0014±0.0001 

0.07 

38.12 

41.90 

-3.78 

MRK 766 

1.98±0.17 

56.3 

0.0075±0.0006 

0.16 

39.10 

43.28 

-4.18 

NGC 5548 

1.6±0.3 

75.0 

0 .011±0.002 

0.19 

38.76 

42.93 

-4.17 
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Table 3. Black hole mass, size of radio-sphere in gravitational radii, and bolometric luminosity 


Object 

T/bh 
10® Mq 

Ref (Mbh) 

Rpc 

log -Lbol 
ergs“^ 

Ref (Lboi) 

Lbol/LEdd 

MR2251-178 

200 

Lira et al. f2011j 

38 

45.3 

Lira et al. f2011j 

0.08 

NGC 3783 

29.8 

Peterson et al. ("20041 

52 

44.4 

Woo & Urrv 120021 

0.07 

NGC 5506 

2.5 

Naear et al. (2002) 

710 

44.1 

Soldi et al. ('2011) 

0.4 

NGC 7469 

10 

Peterson et al. ("20141 

330 

45.3 

Woo V Urrv 120021 

1.5 

ARK 564 

2.6 

Botte et al. (2004) 

930 

45.0 

Romano et al. (2004) 

3.0 

NGC 3227 

42.2 

Peterson et al. 120041 

17 

43.9 

Woo V Urrv 120021 

0.014 

MRK 766 

1.8 

Bentz et al. (2009) 

920 

44.0 

Bentz et al. (2009) 

0.4 

NGC 5548 

67.1 

Peterson et al. ("20041 

30 

44.8 

Woo & Urrv 120021 

0.08 


Mbh is taken d irectly from the lite rature, except for NGC 5506 in which it is estimated by using the broad Paschen line width and 
luminosity from lNagar et ^ ||2002| V 


the high-resolution A-array configuration. Table[4]lists these 
archival core flux density measurements, as well as the cur¬ 
rent 95 GHz measurement. Since the current measurements 
have larger beams than those of the VLA A-array in the 
different bands, we also present in Table [5] archival flux den¬ 
sity measurements obtained with other VLA configurations 
which basically have larger beams than in Table |4l 

The derived radio broadband Spectral Energy Distri¬ 
butions (SEDs) are plotted separately for each source in 
Figure [2l including also the present detection at 95 GHz. 
One must bear in mind that the unmatched beam sizes at 
different frequencies may lead to inconsistent radio spectra 
if different scales of radio emission are present, such as an 
extended component from a putative jet. However, large flux 
discrepancies between different beams that have an appre¬ 
ciable effect on the spectrum are found only for MR 2251-178 
and for NGC 5548 (see Figure [2)), indicating the flux of the 
other sources likely originates in the inner sub-arcsec region 
of the AGN. Moreover, since the data in Table |4] are taken 
over many years, and these sources have been observed to 
vary, the SEDs in Figure [2] should be treated with caution. 
Below we describe the observed spectral pattern found for 
these sources. 

In Fig. H we show the spectral slope indices a (Eh oc 
zz““), obtained from the low-frequency flux densities, ex¬ 
cluding the 95 GHz data point. When available, we provide 
one slope for the VLA A-array flux density (circles), and one 
for the larger beams (triangles). The spectra for the most 
part are steep with 0.8 < a < 1.3. Mrk766 has a somewhat 
flatter spectrum with a = 0.68 when the large beam re¬ 
sults are considered. MR 2251-178 and NGC 5548 have flat¬ 
ter spectra (a — 0.39,0.14, respectively), when the NVSS 
(45”) flux density at 1.4 GHz is ignored. For the most part, 
the 95 GHz flux density overshoots the low-frequency slope. 
The excess factor for each source is given at the rightmost 
column of Table IH providing a range when two slopes with 
different beam sizes are available. The excess is large (x4—7) 
for NGC 3783, ARK 564, and NGC 3227, and smaller (but 
still > 1) for MR2251-178, NGC 5506, and NGC 7469. For 
Mrk766 and NGC 5548 the picture is murkier, although 
some excess can not be ruled out. Apart from the mere de¬ 
tection, this 95 GHz excess is the main result of the present 
campaign. 

It is difficult from our data to pinpoint the frequency at 
which the spectrum turns over from steep to fiat or inverted. 


lAntonucci fc Barvainid (Il988l 'l and iBarvainis et al.l (Il996l ~l 
identified spectra that turn over already at 10- 15 GHz. 
They too ascribed this to a self absorbed synchrotron com¬ 
ponent, and termed it the High-Frequency Excess. In the 
present sample, the sources that have archival data above 
10 GHz are only NGC 5506, NGC 7469, and NGC 3227. Only 
NGC 7469 shows a possible turnover below the 95 GHz data 
point (see Fig. [2]). In contrast with those previous works, in 
which some sources had a high-frequency excess and some 
did not, the present sources may all have at least some 
95 GHz excess with respect to the lower-frequency slope. 


3.3 Beam sizes and variability 

There are caveats to these results, however. Since VLA beam 
sizes comparable to those at 95 GHz (2”2 -5”5 here) are 
not available at all frequencies, and in particular since the 
15 GHz and 22 GHz observations are at sub-arcsec resolu¬ 
tion, the magnitude of the 95 GHz excess reported above 
would change if there is extended emission on scales of a 
few arcseconds, but not if the 95 GHz flux is genuine core 
emission, as suggested by the radio-sphere sizes in Table [S] 
Moreover, since the observations are often many years apart, 
the 95 GHz excess needs to be confirmed with simultaneous 
multi-frequency observations. 


Radio var i ability of RQ AGN has been well documented 


(IWrobell 2000l: Barvainis et al.l 20051; 

Anderson & UlvestadI 

I 2 OO 5 I: Mundell et al.l 20091: Jones et al 

l201lll.lWrobe]l(200nll 


monitored NGC 5548 at 5 and 8.5 GHz and found variability 
of up to ~52% on time scales of months to years, which is 
reflected in the larger error bars for this source in Fig. [2] 
and which is much larger than individual measurement un¬ 
certainties (plotted for the other sources). We further note 
that NGC 5548 is in a historic ally low X-ray and UV flux 
state for the past few years dKaastra et al.l |2014| L If the 
95 GHz emission is coming from similar inner regions as 
the X-ray and UV, it is not surprising that the 95 GHz flux 
now is much lower than expected. Given that all of the low- 
frequency measurements of NGC 5548 are from the 1990s 
(see references in Table [4|, we suspect NGC 5548 would fea¬ 
ture a stronger 95 GHz excess if all its radio measurements 
were simultaneous. In summary, it would be useful to con¬ 
firm the 95 GHz excess for the entire present sample with a 
contemporary multi-frequency campaign. 
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Figure 2. Broad band radio spectral distributions of the eight detected AGNs. 95 GHz flux densities (squares) are from the present 
survey, while lower frequency measurements (circles) are taken from the VLA, as much as available. Note that configuration and beam 
size vary dramatically between observations. The dashed line is the slope (value indicated) fitted to the high-resolution A-array flux 
densities (circles). Where available, a dashed-dotted slope is provided for flux densities obtained with larger beam size (x). In both 
cases, the 95 GHz measurement is excluded from the fit in order to give an idea of the 95 GHz excess. Error bars represent measurement 
uncertainties, except for NGG 5548 where they reflect variability. 
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Table 4. Array-A VLA flux densities in mjy. Excess factor at 95 GHz is relative to the low-frequency spectrum. 


Frequency 

Resolution 

1.4 GHz 

0"3 

5 GHz 

0'.'33 

8.5 GHz 

O'.'20 

15 GHz 

0'.'13 

22 GHz 

0"089 

95 GHz 

2:'2-5'.'5 

Object 


Ref 

F. 

Ref 

F. 

Ref 

F^ 

Ref 

Fu 

Ref 

F^ 

Excess 

MR2251-178 

S.OiO.l 

map 

S.liO.l 

map 

— 

— 

— 

— 

— 

— 

1.6±0.3 

1.7-27.7 

NGC 3783 

43.6±2.0 

NVSS 

13.0±1.0 

UL84 

7.7±0.2 

SCOl 

— 

— 

— 

— 

3.1±0.7 

4.1 

NGC 5506 

304±5.0 

ORIO 

136±4.0 

ORIO 

84±4 

ORIO 

44±5.0 

UL84 

22±4 

TAll 

10±1 

1.2-1.3 

NGC 7469 

134±6.0 

UN87 

26.9±2.0 

LA04 

14.8±0.3 

ORIO 

10.6±0.3 

ORIO 

17.5±0.5 

PR04 

5.0±0.2 

1.3-1.7 

ARK 564 

25.7±0.2 

map 

8.6±0.4 

LAll 

3.0±0.5 

SCOl 

— 

— 

— 

— 

1.1±0.2 

4.8 

NGC 3227 

78.2±5.0 

HOOl 

15±2 

UL81 

12.2±1.3 

KU95 

4.7±0.3 

NA05 

— 

— 

4.1±0.2 

5.3-6.8 

MRK 766 

36.4±1.9 

UL89 

11.7±0.6 

UL95 

8.7±0.4 

KU95 

— 

— 

— 

— 

2.0±0.2 

1.0-1.8 

NGC 5548 

5.1±0.5 

HOOl 

4.2±0.4 

WROO 

3.9±1.8 

WROO 

— 

— 

— 

— 

1.6±0.3 

0.6-2.8 


Refere nce (b eam size indica t ed in parenth esis when deviates from nomina l): NVSS (45^ ^ l[Condon et al.llll998h. UL84lUlvestad Sz. WilsonI 
lll9^), SCOllSchmitt e t alj ll2 00ll), OR IO lOrienti Sz Priet^ ||2010| ), T All iTarchi e t alj <l20 1l[). UN 87 lUnger et al ] ||l987|) , LA04 (1") 
iLal et ^ l|200 ^. LAll | Lal et alJ II2011]). ULSl jlRve^ ad et al.N l98l|) . KU95 | K ukul^ et al.l l|l995h . P R04 (2(^5 ) IProuton et alJ l|2Q0‘A . 
NAOs iNagar et alj ||2005h . HOOI IHo &: UlvestadI ||2001^ . UL89 lUlvestad &; WilsonI lll989h and WROO [Wrobell ||2000|) . ’ map’ stands for 
the map obtained from the VLA archive; For MR 2251-178 map beam sizes are 676 at 1.4 GHz and l78 at 5 GHz. 


Table 5. Flux densities in mJy and beam size in arcseconds from VLA observations with larger beam sizes than in Table [4] 


Frequency: 

1.4GHz 

5 GHz 

8.5 GHz 

15 GHz 

Object 

F^ 

Ref (beam) 


Ref (beam) 


Ref (beam) 


Ref (beam) 

MR2251-178 

16.2±0.1 

NVSS (45") 

— 

— 

— 

— 

— 

— 

NGC 3783 

— 

— 

— 

— 

— 

— 

— 

— 

NGC 5506 

337.9±0.1 

FIRST (5") 

139.2±0.1 

C096 (3") 

— 

— 

58±3.0 

ORIO (077) 

NGC 7469 

145.7±0.1 

FIRST (5") 

66.1±3.3 

ED87 (15") 

42±4.0 

KU95 (275) 

— 

— 

ARK 564 

28.6±0.2 

NVSS (45") 

10.7±0.2 

LA04 (1") 

— 

— 

— 

— 

NGC 3227 

82.8±0.1 

FIRST (5") 

35.0±0.7 

GA06 (15") 

16.8±0.5 

KU95 (2'.'5) 

4.44±0.6 

map (078) 

MRK 766 

40.3±0.2 

FIRST (5") 

20.4±0.6 

GA06 (15") 

10.3±0.4 

KU95 (2'.'5) 

— 

— 

NGC 5548 

24.4±0.2 

FIRST (5") 

10.5±0.5 

ED87 (15") 

4.4±0.3 

KU95 (2'.'5) 

— 

— 


Reference: NVSS iCondon e t al.l lll 99^ ), FIR ST [Becker et alj l|l99f ). C096 [Colbert_et_ajj _jl99d). ORlo lOrienti & Prietol ll2010l) . ED87 
lEdelsor] <ll987l) , KU95 lKukuh^t al.l 1 19951) , LA04 lLa^^dT ||20^ ), and GA06 iG'alhmo^ et alj <l2006l) . ’ map’ stands for the map 
obtained from the VLA archive. 


4 DISCUSSION 

Since none of onr sonrces is resolved, we need to base our 
analysis of the core emission on the spectrum. On one hand, 
a steep SED with a > 0.5 is typical of optically-thin syn¬ 
chrotron emission, while a flat SED (a < 0.3) is interpreted 
as the superposition of several optically-thick regions, where 
the synchrotron emission is self-absorbed. Steep spectra have 
been customarily associated with extended jet emission on 
a kpc scale, while flat spectra likely come from a pc-scale 
core. These behaviors have been commonl y used to inter¬ 
pret both RL and RQ AGN emission (e.g., iBarvainis et al.l 
llQOdf . although the power and speed of the jets, if they can 
be cal led jets, in RQ AGN is believed to be su b -relat ivistic 
(e.g., iFalcke et al.l Il99fil: Ulvestad et all llOOSl . l2005h . Re- 
centlvT lKharb et al.l ||2015 1 used parsec-scale VLBI resolu¬ 
tion to suggest that these steep-spectrum jets in RQ AGN 
are in fact coronal winds. The low-frequency radio emis¬ 
sion from RQ AGN is a mixed bag of emission mecha- 
nisms, including complex morphology and varying spectra 
([Baiwajnis_et_alJ[ljl9fi ^^ et al’lflQ^ : iGallimore et al.l 

l2006l : Giroletti fc Panessa 2009h . H^e. we discuss the pos¬ 


sibility, based on the small Rpc sizes that the 95 GHz cores 
are due to the accretion disk or to outflows in its vicinity 
that are not jetted in nature, namely are neither relativistic 
nor well collimated. These cou ld be similar to stel l ar coro¬ 
nal m ass ejections (CMEs; e.g.. lBastian et ^Il998l : iBastianI 
l200itl . but on the larger scales of an AGN accretion disk. 

A similar effect of 95 GHz excess emission (diminishing 
steepness with frequency, an d in some cases even an inverted 
spectrum) was o bserved bv IPoi et al.l (|2005l . 1201 il l and by 
IPark et al.l ll2013ll in lo w luminos i ty AG N as well as in early 
type galaxies. Indeed, IPoi et akl ( ioij) conclude that this 
is indicative of core AGN emission. Extended stellar dust 
emission is also discussed, but is unlikely, as dust emission 
from stars, let alone the AGN, is hotter and becomes impor- 
tant only at hig her frequencies above 300 GHz. Moreover, 
IPoi et al] (I 2 OIIII report variability of several sources on short 
time scales, which rules out a stellar component. Interest- 
ingly, the most c onspicuous 95 GHz excess in the sample of 
IPoi et aP (I2OIIII is found in sources that also reveal very 
broad optical emission lines when observed at high angular 
resolution, namely in NGC2782, NGC4143, and NGC4579 
dShields et al.llioOTl') . 
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The mm-wave emission is hardly absorbed once it leaves 
the radio core, except perhaps by a Compton thick column 
of ionized gas, and since free-free absorption strongly de¬ 
creases with frequency. Consequently, we expect that the 
95 GHz excess will be observed in type II AGN as well. In 
fact, mm-wave surveys could ex ploit the excess to effici ently 
detect obscured AGN. Indeed, iGallimore et al.l ll2004l ') find 
an inverted spectrum for the radio spot (SI) that is believed 
to be the location of the nucleus of the archetypical type II 
AGN, NGG 1068. 

In a n attempt t o exp lain this spectral change around 
100 GHz, iDoi et al] (l201lh invoke advection dominated ac¬ 
cretion flows (ADAF), since similar behavior is observed in 
Sgr A*. However, the present X - ray se lected sample is quite 
different from that of iDoi et al.l ll201lll . who chose their tar¬ 
gets based on their low luminosity. The present sample in¬ 
cludes high Lboi/TBdd sources, as listed in Table [3] (exceed¬ 
ing unity for NGC 7469 and Ark 564). Therefore, their accre - 
tion disks are highly radiative dShakura fc Sunvaevlll973l . 
and ADAF is irrelevant. 

We thus suggest the observed high-frequency excess 
emission is of coronal nature, in the sense that magnetic 
activity around the accretion disk is responsible for the mm 
wavelength component. It originates further in the AGN 
core where the energy density of the electrons is higher 
and the magnetic field is substantially stronger compared 
to the more extended regions, from where the low frequency 
optically-thln synchrotron component arises. The magnetic 
activity of the disk may increase with accretion (Tboi/TBdd), 
and results in increased 95 GHz emission. In the next sub¬ 
section, we further examine the place of our sample in the 
context of other 95 GHz measurements from the literature. 


4.1 Radio to X-Ray Luminosity Ratio 


In order to further study the origin of the 95 GHz emission of 
RQ AGN and their possible connection with the accretion 
disk, we wish to compare with the X-ray emission, which 
is known to arise from the center of the AGN and to be 
unbeamed. For example, if the 95 GHz emission is linked 
with the disk corona, it could manifest itself in a different 
X-ray/radio connection than that of RL AGN, where a rel¬ 
ativistic jet plays a dominant role in the radio emission. 

Starting with RQ AGN, Figure |3] (left panel) presents 
the comparison between the 95-GHz luminosity and the 
2-10 keV X-ray luminosi t y for the present sample, the 
RQ sources of iDoi et all |201lll , and the upper limits 
from lOwen fc Puschelll 1 1982 ), which includes muc h brighter 
source s. We determined the RQ sources from iDoi et aP 
ll201lll based on their core 5 GHz l uminosity and the a ssoci- 
ated logLi{/Lx < —2.8 criterion llPanessa et akllioOTl fl. In 
the figure, we plot only 95 GHz flux densities, and average 
over multi-epoch data when present. More specifically, for 


NGC 4143 the plotted mean flux density is 7.6 m jy, and for 
NGC 4258 it is 8.6 mJy. Many of the RQ targets iniDoi et all 
ll201lll and all of those in I Owen fc Puschelll (Il982l ') are upper 

limits. Detections are plotted as full symbols. _ 

_ It can be seen that the detected sources in iDoi et al] 

ll201lll are concentrated below log La ^ 39 {La = I'Ll, in 
ergs“^ at iz = 95 GHz). This is similar or only slightly 
less luminous than the present sources, except for MR 2251- 
178 which is more luminous (logLij = 40.18). Given that 
log La ^ 39 is about the sensitivity limit of small millime¬ 
ter arrays observing nearb y galaxies, and that many of the 
sources in iDoi et al. (1201 il l have only upper limits, it is un¬ 
clear how significant the differenc e at 95 GHz between the 
samples actually is. The survey of lOwen fc Puschdl (Il982l l 
is much less sensitive, and its RQ sources all have only 
upper limits at log La < 42, so also consistent with all 
RQ AGN lying below log La ^ 40. On the other hand, 
the sources of the present samp le in X-ra y s lie between 
42 < log Lx < 44, while those of iDoi et all (I 2 OIIII occupy 
the log Lx < 42 region o f the s catter plot in Fig. [3l and 
those of lOwen fc Puschelll (|l982h are all X-ray bright with 
log Lx > 44. 

The eight present sources lie within a factor of 2 from 
the La = 10~^Lx relation. Thi s is reminiscent of th e La ~ 
10~®Lx relationshi p at 5 GHz dLaor fc Beharir2008l ') that is 
also found in stars iGiidel fc Ber^ 1993 ). Bear in m ind that 
the X-ray luminosity used bv iLaor fc Beharl (|2008l 'l was in¬ 
tegrated over the 0.2-20 keV band, while in Fig. [3] we take 
the 2-10 keV luminosities, as these are available in the lit¬ 
erature. Assuming a mean X-ray energy spectral slope of 
-1, the present Lx values would increase by a factor of 
~2.9 for the entire 0.2 - 20 keV band. The frequency used 
here is 95 GHz, but for a radio spectral index of a = 1, 
La = vLv is independent of frequency. The relation adjusted 
for the extended X-ray band would therefore yield on aver¬ 
age La 3x 10~®Lx. Furthermore, X-ray absorbed quasar s 
were excluded from the sample of iLaor fc Beharl I 2 OO 8 II . 
while most of the present AGNs are absorbed. Accounting 
for X-ray absorption would bring the present La — Lx rela¬ 
tion even closer to the canonical La 10~®Lx value, which 
holds at 5 GHz. 

The sources from the other samples in the left hand side 
of Fig. [ 3 ] do not appear to lie along the La = 10“‘*Lx re¬ 
lation, although many are upper limits. The eight sources 
in the present sample were chosen based on their X-ray 
brightness and variability. Hence, it is dangerous to reach 
a conclusion based only on their La/Lx behavior and how 
it is different from other samples with other selection crite¬ 
ria. Although their position on the La/Lx diagram (Fig. [3} 
appears unique {La/Lx « 10”'*), this could easily be the 
re sult of the select ion bias. Indeed, sources from the sample 
of iDoi et al. I dioTIl) have lower X-ray luminosities, but radio 
luminosities comparable to the present sample. On the other 
hand, the many upper limits impede a conclusive argument 
about the general La/Lx behaviour. 


^ According to this criterion, the RQ AGN are NGC 404, 
NGC 2768, NGC 2974, NGC 3065, NGC 3156, NGC 3610, 

NGC 4233, NGC 4494, NGC 4742, NGC 5173, NGC 5283, 

NGC 5812, NGC 266, NGC 3147, NGC 3169, NGC 3226, 

NGC 3718, NGC 4203, NGC 4143, N GC 4258, NGC 4565, 

NGC 4579, and NGC 4772. The rest of the lDoi et al.l ll201lh sam¬ 
ple are RL AGN. 


4.1.1 Radio Loud Samples at 95 GHz 

Expanding the picture to include RL AGN, we plot in the 
right hand panel of Figure |31 95 GHz and X-ray measure¬ 
ments of FRI and FR II low and h i gh ex ci tation galax- 
ies (LEG and HEG) from iDoi et ahl I2OI ih . ICotton et al.l 
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Figure 3. Lij(= p'Li/) at 95 GHz vs. Lx (2—10 keV) in erg comparing the present sample (diamonds) with published observations of 
RQ {left an d right') and RL ( right) AGN cores. For RL sources, only detections are plotted. The present sample is more X-ray luminous 
than that of iDoi et al.l (l201lh and tightly follows the Lji = 10“^Lx correlation, while the RL LEG sources are spread around Lji = Lx- 


ll2009l . CO in the figure), and iHardcastle fc LoonevI ll2008l . 
H&L in the figure), who observed a small sample of pow¬ 
erful radio cores at 90 GHz. We also include the bright 
sample (> IJy) of RL quasars wit h flat or inverted spec- 
tra at millimeter wavelengths from I Owen et al.l ([loU) and 
th e mixed sample of RQ an d RL quasars observed at 90 GHz 
bv lOwen fc Puschelll (Il982l) . As anticipated, the RL sources 
occupy the upper envelope of the plotted distribution, which 
corresponds to higher radio luminosities at given Lx values. 

Among the RL AGN, a tight linear correlation is found 
for the LEG galaxies, both FRl and two FRII radio sources. 
Considering the censored data, the relation is consistent 
with Lr/Lx = 1.0 ± 0.1, with a very low probability 
of < 0.0001 that this is fortuitous. This relation com¬ 
plements similar relations found for LEG sources between 
the radio core power (i.e., at 5 GHz) and multi-band nu- 


(e.g.. Chiaberge et al.l 199f 

. l200i Hardcastle & Worralll 

200nl: iBalmaverde & Gaoetti 

I 2 OO 6 I: IBalmaverde et al.ll200fil: 

Baldi & Gaoettil 2010l). These relations are internreted as 


evidence for synchrotron emission from the basis of the jets 
at different frequencies. The Lr/Lx = 10“"^ values of the 
present sample are different by 4 orders of magnitude sug¬ 
gesting its underlying physical mechanism is completely dif¬ 
ferent. Having said that, the basis if the jet and an accretion 
disk corona are not necessarily mutually exclusive. 

On the other hand, the RL HEG galaxies are brighter 
in the X-rays than the LEG by two orders of magnitude or 
more (Fig. |3}. This X-ray emission is interpreted as coming 
from the a ccretion disk, in excess of t he beamed compo¬ 
nent (e.g.. IHardcastle fc Worralll Dood : iDicken et al.l l2009l : 
iBaldi et al.fl201Clh . and is the dominant X-ray source in RQ 
AGN. At yet higher 95 GHz lumin osities, Fig. [3 inclu des also 


the flat-spectrum radio quasars of lOwen et al.l lll98ll) , which 
are even slightly more p owerful in the radio than those of 


are even sligntly more p o 
lOwen fc Puschelll (Il982l) . 


Although RL AGN are generally brighter at 95 GHz 
than the RQ AGN, the right hand panel of Fig. [3] lacks 


a clear bi-modality between RQ (Lr/Lx <C 1) and RL 
(Lr/Lx ~ 1) AGN. Especially at low luminosities the dif¬ 
ference between RL and RQ decreases, and the two sam¬ 
ples appear to approach each other on the lower Lr and 
Lx corner of Fig. [S] This could be a result of the mixed 
and incomplete sample presented in the figure. On the other 
hand, if the black-hole mass is a crucial property f or ra¬ 
dio loudness, and high black - hole mass of >10^~^ Mrx llLaorl 
l200d : iBaldi fc GaDettil|20ld : IChiaberge fc Marconill201ll ) is 
required to launch relativistic jets, then the continuous BH 
mass values in AGN could explain the continuous range of 95 
GHz luminosities in the figure, where the RQ AGN, which 
generally do not exceed 10'*° ergs“^ represent a 95 GHz 
disk contribution with no jets. Large scale (kpc) radio jet 
emission tends to produce steep spectra, contributing most 
significantly at low frequencies. The difference therefore be¬ 
tween RL and RQ is more prominent at low frequencies than 
at 95 GHz. 


5 CONCLUSIONS 

We observed eight X-ray bright and variable RQ AGN at 
95 GHz with GARMA and ATCA. All of the sources were de¬ 
tected. Millimeter-waves are the sole remaining unexplored 
band of AGN physics, and the present observations help to 
fill this band gap. The emerging SED, based on archival data 
shows significant mm-wave emission from the inner core of 
the AGN, which dominates the spectrum at high radio fre¬ 
quencies. 

The main features of the 95 GHz emission discovered in 
the present work can be summarized as follows: 

• The observed sources show 95 GHz excess of up to x7 
with respect to their extrapolated steep spectra at lower 
frequencies. 

• The deduced size Rpe for an optically thick 95 GHz 
radio-sphere is 10“"* — 10“^ pc, and in terms of gravitational 
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radii ranges roughly from 10 - 1000 Tg, similar in size to that 
of the X-ray source. 

• All eight sources lie tightly along an Lr — 10~‘^Lx 
correlation. 

We interpret these findings as accumulating evidence 
for the 95 GHz emission being associated with the accre¬ 
tion disk, more specifically we suspect it originates near the 
X-ray corona. The excess emission implies a compact, op¬ 
tically thick core. The Lr/Lx correlation is analogous to 
that found for RQ AGN and for stellar coronae at 5 GHz. 
These are two separate lines of (circumstantial) evidence 
for an accretion disk/coronal origin of the compact 95 GHz 
source that is physically different from the majority of radio 
emission in RL AGN. 

We suggest following up the present exploratory cam¬ 
paign of an X-ray selected sample with a systematic study of 
the SEDs of RQ AGN between 1 - 100 GHz for a complete 
sample, selected independent of radio and X-ray properties. 
For example, in the PG quasar sample this will allow to 
compare with the vast data available, and to explore possi¬ 
ble dependence of the mm flux on the continuum and line 
emission properties in other bands. 

In order to further investigate the conjecture of radio 
emission from the accretion disk corona, where the X-rays 
originate, simultaneous millimeter and X-ray observations 
are needed. Temporal correlation between the X-ray and 
radio/millimeter light curves from RQ AGN would be the 
smoking gu n of coronal ra dio emission. For example, a Neu- 
pert effect (lNeuDertlll968h , as first observed o n the Sun, and 
later in stellar coronae fe.g.. lGudel et al.l2002|j is regarded as 
direct evidence for chromospheric evaporation that produces 
X-ray flares, as a result of coronal heating by non-thermal 
electrons that produce the radio emission. This effect, or 
similar, in AGN would constitute more evidence that the 
radio emission in RQ AGN is of coronal origin, and linked 
to the X-ray emission from the immediate vicinity of the 
accretion disk. 
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